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Abstract 

The self-consistent separable RPA (random phase approximation) method is formulated for 
Skyrme forces with pairing. The method is based on a general self-consistent procedure for fac¬ 
torization of the two-body interaction. It is relevant for various density- and current-dependent 
functionals. The contributions of the time-even and time-odd Skyrme terms as well as of the 
Coulomb and pairing terms to the residual interaction are taken self-consistently into account. 
Most of the expression have a transparent analytical form, which makes the method convenient 
for the treatment and analysis. The separable character of the residual interaction allows to avoid 
diagonalization of high-rank RPA matrices and thus to minimize the calculation effort. The previ¬ 
ous studies have demonstrated high numerical accuracy and efficiency of the method for spherical 
nuclei. In this contribution, the method is specified for axial nuclei. We provide systematic and 
detailed presentation of formalism and discuss different aspects of the model. 


1 


I. INTRODUCTION 


Effective nuclear forces (Skyrme, Gogny, ...) are widely used for description of diverse 
properties of atomic nuclei (see, for review [lj). However, their application to nuclear dy¬ 
namics is still limited even in the linear regime which is usually treated within the random- 
phase-approximation (RPA). The theory is is plagued by dealing with high-rank matrices 
which make the computations quite expensive. This is especially the case for non-spherical 
systems with their demanding configuration space. The rank of the matrices is determined 
by the size of the one-particlc-one-holc (1 ph) space which becomes really huge for deformed 
and heavy spherical systems. 

RPA problem becomes much simpler if the residual two-body interaction of a given mul¬ 
tipolarity A/i is factorized (reduced to a separable form): 

K 

< mn \ V res\V > a+a+ajdi -> ^ ^ k ,X^X^, 

mnij k,k'=l 

A" = E < pAM* > (i) 

ph 


where X k ^ are hermitian one-body operators and K^ k , are strength constants. The fac¬ 
torization allows to reduce a high-rank RPA matrix matrix to a small matrix with a rank 
determined by the number of the separable terms. The main problem is to accomplish 
the factorization self-consistently, with minimal number of separable terms and with high 


accuracy. 

Several self-consistent schemes U 


were proposed during last decades 


and signified a certain progress in this problem. However, these schemes are not sufficiently 


general. Some of them are limited to analytical or simple numerical estimates 


the others are not fully self-consistent 17]. Quite promising is the approach j8j, [9f] for Skyrme 
forces. However, it still deals with RPA matrices of rather high rank (~ 400). Besides, 
it neglects contributions to the residual forces from the Coulomb interaction and time-odd 
densities and currents . 




In this connection, we proposed some time ago a general self-consistent separab 
(SRPA) approach relevant to arbitrary density- and current-dependent functionals 


EEa]. 


e RPA 


lfl 11 

The method was implemented to the Skyrme functional [lj, |U| . In SRPA the one- 


body operators and strength constants of the separable expansion are unambiguously derived 
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from the given energy functional. Since we use the self-consistent procedure, there is not 
need in any new parameters in addition to those in the initial functional. The factorization 
dramatically reduces rank the of RPA matrix. Usually a few separable terms (or even one 


term) are enough for accurate reproduction of the genuine residual interaction jloj|. Such 
impressive result becomes possible due to the effective self-consistent procedure [4( based on 
solid physical arguments. Besides, the separable operators are constructed to have maxima 


at different slices of the nucleus and thus cover both nuclear surface and interior [Q. Hence 
SRPA exhibits accuracy of most involved RPA versions but at much less expense. 

One of the main SRPA advantages is its simple and transparent formalism which makes 
the method very convenient for the analysis and handling of the numerical results. Being 
self-consistent, SRPA allows to identify spurious admixtures connected with violation of the 
translational or rotational invariance. 

SRPA exploits the full lph space and thus equally well treats collective and non-collective 
states. However, quite often (e.g. for giant resonances) we do not need the detailed RPA 
description. Then the strength function method, when we completely avoid the calculation 
of the RPA states, is much more optimal. The separable character of SRPA allows to 
construct the strength function based on the Lorentz smoothing function. The strength 
function is naturally separated into two terms, from the mean-held and correlations. 

SRPA has a peculiarity to incorporates to the residual interaction the contributions of 
both time-even and time-odd variables (densities and currents). The time-odd variables 
naturally appear in the Skyrme functional constructed to contain all the possible bilinear 
forms from the basic nucleon and spin densities together with their derivatives up to the 


16]. 


second order ll6|. Besides, time-odd densities and currents are necessary to keep the Galilean 
and gauge invariance of the Skyrme functional fiol . Ir3 |. Skyrme functionals of this kind are 
actively used for investigation of both ground state and dynamics of atomic nuclei (see e.g. 
[ul j and references therein). 

The recent studies with Gogny forces show that contributions of the spin-orbital and 
Coulomb forces to the residual interaction can be important for the description of low- 
lying states and giant resonances in exotic nuclei ju|. SRPA takes into account both these 
contributions. We plan to scrutinize their impact in our further studies. 

SRPA is quite general and in principle can be applied to a variety of finite Fermi systems 


and different functionals. For example, it was derived for the Kohn-Sham functional 


W2Q] 
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and widely used for descri 
deformed atomic clusters 121 


ion of linear dynamics of valence electrons in spherical and 


22 


m 


m 


m 


In the previous SRPA studies, we considered nuclear dynamics in spherical nuclei QQ. 
However, ability of SRPA to reduce the computational effort is much more decisive for 
deformed systems with its huge lph configuration space. So, in the present paper we specify 
the Skyrme SRPA formalism to axial atomic nuclei. The pairing and its contribution to the 
residual interaction are taken into account. One of the aims of this paper is to present the 
SRPA by systematic and even tutorial way, with all necessary details. 

The paper is organized as follows. In Section |Hl the derivation of the general SRPA 
formalism is given and discussed. The method is SRPA is specified for Skyrme forces in Sec. 
11111 The choice of the initial operators is discussed in Sec. IIVI The summary is done in Sec. 
El Details of the formalism can be found in Appendices A-H. 


II. BASIC SRPA EQUATIONS 


A. Main requirements 


The present model provides the self-consistent factorization of the residual interaction to 
the explicit form 

1 K 

^ Koo “ V V \ I^Rk.R'k 1 X R kX^'k 1 ~\~ Tjsk s'kiYskYs'k 1 f • (2) 


2 ^ ^ ^ {^skjS'k'X s kX s 'k' T ?7 sk,s'k'Y s k^ s'k'} ■ 

ss' k,k '=1 


We assume the residual interaction of a fixed multipolarity (A for spherical nuclei and Ap for 
deformed nuclei) but, for simplicity, skip hereafter the multipole index. In (( 21 ), indices s and 
s' label neutrons and protons; X sk and Y sk are time-even and time-odd hermitian one-body 
operators. Their time-parity properties formally read 

TX sk T 1 = l s k, 7 * — +1, 

TY sk T- 1 = 7 %Y sk , It = -1, 

where T is the operator of time inversion. The expansion © includes time-odd operators 
oecause some Skyrme functionals possess both time-even and time-odd densities and currents 
11, see their list in the Appendix O As was mentioned above, the time-odd densities are 

nn 

necessary to keep the Galilean and gauge invariance of the Skyrme functional [15|, 1161] . 
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Though these densities do not contribute to the static mean field Hamiltonian of spin- 
saturated systems, they can provide time-dependent perturbations and thus have to be 
taken into account in the description of nuclear dynamics. 

The presence of both time-even and time-odd variables naturally leads to formulation of 
the model in terms of hermitian operators with given time-parity. These operators have the 
useful property that 

<0|[i,H]|0>~(l- 7 ^7x) (3) 

i.e. the average commutator does not vanish only for operators A and B having the opposite 
T-parities ( 7 ^ = — 7 ^), see Appendix IP> 21 for more details. This property will be widely 
used in our derivation. 

The model should satisfy some principle requirements. The expansion © has to be self- 
consistent. It should involve the minimal number of the separable terms and, at the same 
time, accurately reproduce the true residual interaction. The operators Xk and Y* and their 
weights should have simple and physically transparent structure. Below we will develop the 
scheme which fulfills these requirements. 


B. Time-dependent Hamiltonian 


The nucleus is assumed to undergo small-amplitude harmonic vibrations around Hartree- 
Fock (HF) or Hartree-Fock-Bogoliubov (HFB) ground state. We start with a general time- 
dependent functional E(Jf(r,t)) involving a set of arbitrary neutron and proton densities 
and currents Jf{f,t) (where s — n,p ; a labels densities and currents) 





(4) 


where |T(f) > is the wave function of the vibrating system described as the time-dependent 
Slater determinant. Time-dependent densities and currents are determined through the 
corresponding operators as 

occ 

J?(f,t) =S< '&(t)\j?{r)\'S(t) > s = J2‘Ph(,r,t)j a (r)(p h (r,t) (5) 

hes 

where <p h (r,t) is wave function of the hole (occupied) single-particle state. The set © 
includes both time-even and time-odd densities and currents, see examples in the Appendix 

□ 
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In the linear regime, the time-dependent density reads as a sum of the static part and 
the small time-dependent perturbation: 




( 6 ) 


Then, substituting © into 0 and keeping the terms up to the linear order for <5J“(r, t), 
one gets the single-particle Hamiltonian 


with the static mean-field part 


h{t) T h res {t) 




(7) 


8J? 


( 8 ) 


and the time-dependent response 


KM = = - 

a.'s' s ols a's' s s 


(9) 


The later determines oscillations of the system. For the brevity of notation, we skip the 
dependence on space coordinates in and hereafter in this section. The explicit space 

dependence of the key SRPA operators and values can be found in the Sec. Hill 

The next step in our derivation is to specify the unknown density variations (or transition 
densities) 

sjfit) =< > - < oi j s q io > (io) 

where |0 > is the static ground state. For this aim we should define the perturbed many-body 
wave function |T(f) >. 

C. Scaling perturbed wave function 


The macroscopic perturbed many-body wave function |T(t) > s is obtained from the static 
HF or HFB ground state |0 > s by the scaling transformation 4] 


K 


!*<*».=n exp[-iq sk (t)P sk \exp[-ip sk (t)Q sk } |0 


>< 


( 11 ) 


k =1 
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Here both |T(t)> s |0 > s are Slater determinants; Q s k(f) and P s k(r) are generalized coordi¬ 
nate (time-even) and momentum (time-odd) hermitian operators with the properties 


Qsk — Q~lki 7 t — 
Psk i Qsk\ph Pski 7 t ^ 


( 12 ) 


where H = h 0 + U res is the full Hamiltonian embracing both one-body and two-body parts. 
The subscript ph in the commutator means the mapping into particle-hole domain. If the 
functional includes only time-even densities, then V res does not contribute to the commutator 
and so one may use ho instead of H. 

Operators m generate time-even and time-odd real collective deformations q s k{t) and 
Psk(t). Using (ED, the transition densities m are expressed through these deformations as 


=< *(t)|.7f|*(t) > - < 0| Jf |0 >= 

- iyifot(t) < 0|[h», Uio > +Psk(t) < 0|[Qst, J“]|0 >} 


(13) 


Then the response Hamiltonian Q can be recast as 


^res (t') ^ T Pski^^^sk} 


sk 


where all time-independent terms are collected in the hermitian one-body operators 




5 2 E 


s' s' a'a a ° 

*- = E* = E‘E \j§r r 

s' a'a b * 


< 0|[At,/f]|0 > E, 


j.j<o|[Q,»,y]|o> J$ 


with the properties 


(14) 


(15) 

(16) 


X, t = A+, = +1, A' 1 = A, 


Y, t = Y,l, = - 1 , A* = -y. 


(17) 

(18) 


As is shown below, X s k and Y s k are just the time-even and time-odd operators to be exploited 
in the separable expansion 0- Following the property 0, time-even densities contribute 
only to X s k while time-odd densities only to Y s k- The upper index s' in the operators ed- 
ED determines the isospin (proton or neutron) subspace where these operators act. This is 
the domain of the density operator Jf! entering ED-ED- 
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To complete the construction of the separable expansion m, we should still determine the 
matrices of the strength constants n s k,s'k' and rj sktS ’ k >. This can be done through variations 
of the basic operators ©: 

SX sk (t ) =< V(t)\X sk \y(t) >-< 0|X sfc |0 >= (19) 

= ^ ^ ] Qs'k'jt) < 0| [Ps'k'i X sk ] |0 >= — y ^Qs'k'{t)K s , k ,_ sk , 

s'k' s'k 1 


SY sk (t ) =< V(t)\Y s k\V(t) >-< 0|F sfc |0 >= 

= 7 ^ , Ps’k'(t) < 0 | \Qs'k'i X sk ] |0 >= — ^ ^ Ps'k'if) r ] s i k i_ s k 

s'k' s'k' 


( 20 ) 


where 


^ s'k'sk 


— K sk,s'k’ ~ 7 < > — 

= /' °ii^. >< h'lio >, 

mm 1 s $ 


( 21 ) 


^s'k' ,sk 


— 7 lsk,s'k' ~ 7 < 0 | [Qs'k',Y s 

r 5 2 E 


/fe] 


> 


( 22 ) 


— 1 dr 




<o\[Qsk,J?]\oxo\[Qs'k>,J$ 


> 


Eqs. (EID-(E2D represent elements of the symmetric matrix which is inverse to the matrix of 
the strength constants in ©• Indeed Eqs. (0-0 can be recast to 


^ ^ s'k',sk$Xsk{t) Qs'k'{t ) ? 

sk 

(23) 

^ 'Vs'k',sk5Ysk(t) Ps'k'i[t) • 

sk 

(24) 


Then the response Hamiltonian (ED gains the form 

^res (t) ^ ^ ^ s'k',sk$X s k{t)X s f k' “I - ^s'k',sk$Y s k{t)Y s f k'} (25) 

s'k' sk 

which leads to the same eigenvalue problem as the separable Hamiltonian 

H = h 0 + V r fv 

with from ©. See also j)|) for relevant discussion. 


( 26 ) 




In principle, we have already in our disposal the formalism for linear regime of the col¬ 
lective motion in terms of collective harmonic variables 


q s k(t) = q sk cos (cut) = ^q s k(e lujt + e wt ) , (27) 

Pak(t) = p sk sin{ujt) = Psk(e luJt - e ~ luJt ) . (28) 

Indeed, Eqs. ED, (HE!), (ED, and ED deliver the one-body operators and strength matrices 
which we need for the separable expansion of the two-body interaction. The substitution of 
the response Hamiltonian ED and the perturbed wave function ED into time-dependent 
HF equation 


4w) >= (ho + k»{t))Mt)> 

at 


(29) 


would result in the eigenvalue problem. The number K of the collective variables (and thus 
of the separable terms) depends on the accuracy we need in the description of collective 
modes (see discussion in Sec. IIVI) . For K = 1, the method reduces to the sum rule approach 

n 

with one collective mode [22! • F° r ^ > 1, we have a system of K coupled oscillators and 


the method is reduced to so-called local RPA |22D |28l suitable for a rough description of 
main branching and gross-structure properties of collective modes. However, the method is 
still not ready to describe the Landau fragmentation. For this aim, we should consider the 
detailed lph space. This will be done in the next subsection. 


D. Coupling with lph space 

Collective modes can be viewed as superpositions of lph configurations. To derive this 
relation, it is convenient to introduce an alternative perturbed many-body wave function 

\m >= (i+E £ i°> < 3 °) 

s phes 

where 

A J = ajah (31) 

is operator of the creation of lph pair and 

<4(i) = cj+e”‘ + (32) 

are time-dependent amplitudes of particle-hole configurations in the perturbed state. Here 
we used the Thouless theorem which establishes the connection between two arbitrary 
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Slater determinants. The wave function 4HD is obviously the microscopic counterpart of the 
macroscopic wave function (HU). 

Substituting © and © into the time-dependent HF equation m, one gets, in the 
linear approximation, the relation between c± h and collective deformations q s ^ and p s k 

S± _ 1 Es'fc'fe'fc' < ph\X s s , k , |0 > Tips'* < ph\Y s s , k ,\0 >] ^oo\ 

Cphes ~ 2 e ph ±u ’ ^ 

where £ p h is the energy of lph pair. 

In addition to Eqs. (USD-©, the variations 5X s }~(t) and 8Y s h(t ) can be now obtained 
with the microscopic perturbed wave function (EDI: 

sx„(t) = J2 E(4w <f*biio> +<„(«) <o\x:’ k \ P h», (34) 

s' phes' 

6V sk {t) = EE «„(*)' <Ph\Yfk\0> +<„(«) <a\Yi\ph>) . (35) 

s' phes' 


E. Eigenvalue problem 


But both amplitudes and collective variables q s k and p s k are still unknown. The time- 
dependent HF equation was already exploited and cannot be used once more to determine 
these unknowns. Thus we need for this aim some additional physical constraint. It can be 
naturally formulated as equality of the dynamical variations of the basic operators 5X and 
5Yk, obtained with the macroscopic m and microscopic (EDD perturbed wave functions. 
Thus we should equate (TBUl - fl^Ul) and (©)-(©)• This gives 


-J2 qs ' k, ^ K s'l’,sk = y y ( c iw* < p h \x s sk\° > + cS ph( t ) <o\x:' k \ P h>), (36) 

s'k' s' phes' 

- &'*'(*) Vs'l',sk = <P h \Ysk\°> +° s ph( t ) < 0|y/fc | ph >). (37) 

s'k' s' phes' 

Substituting pu and (|S2D into these expressions and collecting, for example, the terms with 

e iwt, we g e ^ s 


^ Qs'k'Ks'l^sk 

s'k' 

i P s 'k'Vs'k 1 ,sk 
s'k’ 


2 E E«4T <ph\x:l\o>+c^ <o\x;' t \ P h>}, 

s' phes' 

2 E E {(<*r) - +4 + <o\Y,i\ph> } . 

s' phes' 


(38) 

(39) 
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Then, by using (EH, all the unknowns in the expressions are reduced to the collective 
variables % and pk and we finally get 

V in-\F {XX) - k ~ 1 1 4 .*}-F [XY) \ = n 

/ . X'iskl r s 'k',sk ^sk,s'k'i s'k’,ski U ’ 

sk 


V in -F {YX) + ri ~\F {YY) 

/ s'k',sk ' ^ sk \- s'k’,sk 


^sk,s'k ']} ® 


sk 


(40) 


with 


b43 =EE3-^-^{ < pA|v; t io>*<pAiy%,|o>( £pi +a,) 

s phes Ph 

+ <ph\Xt- k 10 >< 01 X s s , k , |ph > ( £ ph - u )} , (41) 

hS = -*E E ttX q< p h ihtio >*< paivvio > v P „+^) 

s phes £ P h-“ 

+ < ph\Y&\0 >< 0\X s s , k ,\ph > ( e ph - u)}, (42) 

pt fJl = ; E E 5(< phKi 10 >•<PMUVIO > (%* + <V 

, phes £ Ph-“ 

+ < ph|^|/r >< OinVb/r > (e ph -u)}, (43) 

F %Jlk = J2 J2 f 2 * 2 {< ^1^11° >*< P h \Yjk> 1° > 

s phes Ph 

+ < PftlUilO >< 0|p; t ,|pft > (£ pi - w)}. (44) 


These equations can be simplihed using the relations 

< a\Xt- k \ph >=< pftl.YfjlO >•=< pft|.Yf E |0 >= < pft|.Yf E |0 >, (45) 

< 0 Vflpp >=< pp Vi I) >*= - < pp Vi II >= -i< pft|>l|0 > (46) 


which directly follow from the properties (ED and m and expressions (USD and m - Then 
elements of the RPA matrix can be rewritten in terms of real (overline) matrix elements as 


F {XX)_ 
s'k' ,sk 

f (YX) f (xy)_ 

s'k',sk s'k',sk 

F (YY)_ 
s'k' ,sk 


?EE 

s phes 


£ ph < ph\X s .- k \0 > < ph\X s s , k ,\0 > 


£ ph a;2 


-?EE 

s phes 


ph|y||0 > < ph\x s s , k ,\0 > 


£ lh a;2 


^EE 

s phes 


£ ph < ph\YM0 > < ph|y/ fe ,|0 > 


£ ph a;2 


(47) 

(48) 

(49) 


Supposing determinant of the system (jH to be zero, we obtain the dispersion equation 
for determination of the RPA eigenvalues u v . 
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F. Normalization condition 


The standard RPA operators of excited one-phonon states \u > are defined as 

Qi = \ £ Drt - C A) (so) 

s phes 

and fulfill the normalization and orthogonality conditions 


<[<y CJ ]) = Ky, {[Qt,Qt]) = ([4.,4/]) = o, 


(51) 


where A+ h and c v pk are given by m and (BSD , respectively. In the quasi-boson approximation 
{[Aphi A p , h ,\ = Spp'Shh 1 ), the normalization condition [Qu,Qt] — 1 results in the relation 

£ £«ST) 2 - (cg+) 2 } = 2. (52) 

s phes 

By using BSD, it can be recast in terms of the RPA matrix coefficients BZD-BHD: 


££{(<S~) 2 -(<$ + 

phes 


) 2 } 


££ 

s'k' sk 


I Us'k'Qsk 


duj u 


+ 2 q v S 'k'PVk~ 




du v 


+ P V s'k'Pk- 




du v 


(53) 


} = 2 N v . 


Hence the variables q v sk and jf sk should be renormalized by the factor 1 /a/AS- 


G. Pairing contribution 

The pairing Hamiltonian reads 

Vir = ~Y G *xtXs (54) 

s=n,p 

where 

xt = Y a t a ] Xs = Y a ~i a r ( 55 ) 

jes jes 

Then, after the Bogoliubov transformation from particle to quasiparticle operators, the 
hermitian one-body operators have a general form 

A = "Y, < u|A|0 > (afa,j + a^at) (56) 

ij 

Kj,Ki> 0 

= 2 Y {< ij\A\0 > (uiVj + ^UjVi^af at + iT&i&j) 

ij 

+(uiUj - 7rR/^)(b+bj - )} (57) 
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where i are time-reversed states and the time-inverse factor 7 ^ dehnes time-parity of the 
operator A (see details in Appendix IB II) . 

ffence the time-even ( 7 ^ = 1) and time-odd ( 7 ^ = — 1) operators read 

Kj,Ki> 0 

A = 2 < ij|A|0 > (■UiVj + UjVi)(afat + 6qa.j) , (58) 

ij 

Kj,K{> 0 

B — 2 < ij\B\Q > ( u^j — UjVi)(a+at — ajap , (59) 

ij 

i.e. obtain the pairing factors 

u[P = UiVj + UjVi, u[j ) = u^j - UjVi. (60) 

This is the case of time even-operators Q sk and X sk and the time-odd operator Y sk . The 
situation with the time-odd operator 

P sk = i[H, Q sk ] = i{[ho, Q sk } + [l/ r :: p , Q sk ]} = i[h 0 , Q sk ] - Y sk 

is more complicated because of the additional term i[h 0 , Q sk \- Taking into account (|5%j) and 
(??) this operator reads 

Kj,Ki>0 

P sk = 2 {ilwP < W.*|0 > > < ij\YM0 >} (afcA - a&j). (61) 

ijes 

It is seen thut P s k keeps the snme operator structure — ct^CAj us (ED but, at the same 

has, unlike Y sk , the diagonal (i = j matrix elements. The later is the obvious consequence 
of the term i[ho,Q sk \. Explicit expressions for the matrix elements of the operator P sk are 
given in the Appendix IG II 

Besides the factors PI) the pairing results in the specific contribution to the time-even 
part of the response Hamiltonian h res (t). This contribution can be derived in close analogy 
with other time-even densities. Then we get 

ir w = v fc (i) v x‘jr ir) ( 62 ) 

ks s 

where 

X<P atr) = - G s xx,sk ~ v j)( a j a j + a ] a J~)- (63) 

jes 
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The pairing response reads 


Xx,sk = i < Q\[Pks,xt + Xs]|0 > 

= i iejuff < jj\Qks\0 > - vf) < 6|[(a+at - otjotj), (aycp + afaf)} |0 > 

jes Ies 

= 16 YL e i ( u j - v j) M 5i } < jj\Qks\0 > (64) 

jes 

where we used the diagonal part of the operator (HEED- Finally, the pairing contribution to 
the strength matrix is 

[«K]-‘ —i< 6|[n'.,.ht“ ir) ]|6 >= G.xx*Xx ,*• (65) 

It worth noting that in all the sections except of the present one, the pairing factors © 
are supposed to be included into the matrix elements and transition densities and thus are 
not depicted explicitly. 


H. Strength function method 


In exploration of the system response to external fields, we are usually interested in the 
total strength function instead of the responses of particular RPA states. For example, giant 
resonances in heavy nuclei are formed by thousands of RPA states whose contributions in 
any case cannot be distinguished experimentally. In this case, it is reasonable to implement 
the strength function formalism. Besides, the calculation of the strength function is much 
easier. 

For electric external fields of multipolarity EXfi, the strength function can be defined as 


S l (EX^uj) = 

V 

where 

C(w - U V ) = ^ (a; _ a;j/)2 + (A/2) 2 
is Lorentz weight with an averaging parameter A and 

= 4 E <" E < yl/vl° > MT + C T ) 

v s ijes 


( 66 ) 


(67) 


( 68 ) 
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is the matrix element of E\/i transition from the ground state to the RPA state \v >. 
The pairing factor is included to the matrix elements. The forward and backward two- 
quasiparticle amplitudes c^ s± follow from (ITU) (with the subsequent normalization). The 
operator of the electric external held in the long-wave approximation reads 


h, = + Up. (69) 

Further, e e /f is the effective charge (in the dipole case e e J^ = N/A and = —Z/A and e 
is the proton charge); a j u is the energy of the the RPA state \u >. 

It is worth noting that, unlike the standard definition of the strength function with using 
5(w — uj v ), we exploit here the Lorentz weight. It is convenient to simulate smoothing effects. 

For direct use of the expression m , we still have to know all the RPA eigenvalues and 
eigenvectors. This needs an appreciable computational effort. To avoid this, let’s recast (jHHjl 
to the form which does not need the information on the particular RPA states [3D] . For this 
aim, we will use the Cauchy residue theorem. Namely, the strength function will be recast 
as a sum of the residues for the poles z = ±uSince the sum of all the residues (covering 
all the poles) is zero, the residues with z = ±up (whose calculation is time consuming) can 
be replaced by the sum of residues with z = u ± *(A/2) and z = ±£ p h whose calculation is 
much less expensive. Now let’s consider this procedure step by step. 

First, we use (El and rewrite the matrix element (16811 to the form 


Me\ P u 



s ijes 


< ij\f\»\0 > 

r 2 , >2 
£ ij ~ 


^2i £ ijQs'k' < uW'fc'IO > Ai^vV's'w < ufe'10 >} 

s'k' 



V) + frvA^EXn)} 

s'k' 


(70) 


(71) 


where 


A™(E\») = W el" V E J1 < ijjMO > 


A < %(EW=i'£el ,, 'E 


UJ U < ij\Y s s , k ,\0 >< ij\fxu,\0 > 


ijes 


£ ij ^v 


(72) 

(73) 


For the sake of brevity, we introduce the new index f3 = {skg} where g — 1 and 2 for 
time-even and time-odd quantities, respectively. Then the squared matrix element can be 
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written in the compact form 
1 


30 


m 


M W - ^ 


F PP' A A F 00' 


H R 0 R 0' A a A p' = t^Y 1 2N a fi n A 0 A P' = Y1 


00 ' 


2N ^ " -f det D 

00’ 9uj 


aa , det F 

00' du} 


A 0 A 0’ (74) 


where F is the determinant of RPA matrix HDD, Fppt is its algebraic supplement and 


R sk g= 1 — Q s ki R skg=2—Pski (75) 

A,k„,i = A%\E\v), A, k ,, 2 = AZ\E\li). (76) 


Substituting (ITU) to (IQOj) . one gets 


S L {E\^u) 


E 


oo. 


E rpv A v A v 
00’ r 00’ A 0^-0' 

det F u 

OLUv 


t(v 


tO V 


(77) 


Then, taking into account that the determinant det-F^ has first-order poles to = co v , the 
strength function can be rewritten through the residue of these poles on the complex plain, 
or equivalently, through the corresponding contour integrals: 


Sl(E\h\uj) = Res{'A 


,Y,00> F 00'( Z ) A 0( Z ) A 0'( Z ) 

det F(z) 


((uj- z)} z=u 


1 

2 m 



l Y, 00' F 00'(z) A 0(z)Ap'(z) 

Z —-:-—-C(^ 


det F(z) 


z). 


(78) 


Unlike ( 1771 ) . the denominator in (ZED includes the RPA determinant instead of its derivative. 

Following Cauchy theorem, sum of all the residues (covering all possible poles of the 
strength function) is zero and so one can express the residues with z — to u through the rest 
of the others: 


Res[S\ z=U)l/ = -(Res[S] z= -u v + Ses[S] 2 ^oo + Res[S) z=u±i{ A / 2 ) + Res[S\ z=±£ij ) (79) 

where the poles z — to ± i( A/2) and z = Fe i3 originate from the Lorentz weight and 
denominator F PP'( z)Ap(z)Api(z), respectively. The RPA determinant F(z) has zeros only 
at z — ±uj v . 

It’s easy to prove that for L — 0,1, 2 we have lim^i^oo S'(a;, z) = 0. Also, Res[S\ z= - Uv and 
Res[S] z= - £i . can be neglected for large positive values (high energies of giant resonances) 
and relevant values of the averaging parameter A. Remaining residues over the poles z = 
to ± i( A/2) and z = e p h give the final outcome 

Res[S\ z=u , v ~ -Res[S] z=CJ ±i( A / 2 ) ~ R es[S] z=£ij . (80) 
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Finally, the strength function for L = 0,1, 2 reads 


(81) 

z=cj-\-i( A/2) 

+ >2 c(w - £ij). 

s es 

The hrst term in (EH) is contributions of the residual interaction. It vanishes at VJ-es = 0. 
The second term is the unperturbed (purely quasiparticle) strength function. 


Sl(EX/j,lu) = — A 


1 

_ ( 

7r 


z L E,3/j' Fw(z)Ap{z)Ap,{z) 
F{z) 


I. General discussion 


Equations EH-EH, EH-EH, ESI (HOD, EH-©, and © constitute the basic SRPA 
formalism. Before proceeding the specification to Skyrme functional it is worth to comment 
some essential points of the model. 

• One may show (e.g. by using a standard derivation of the matrix RPA) that the separable 
Hamiltonian EH with one-phonon states (©) results in the SRPA equations EH“EH if to 
express unknowns through q k and p k . Familiar RPA equations for unknowns require 
the RPA matrix of a high rank equal to size of the 1 ph basis. The separable approximation 
allows to reformulate the RPA problem in terms of a few unknowns q k and p k (see relation 
EH ) and thus to minimize the computational effort. As is seen from (1401) . the rank of the 
SRPA matrix is equal to 4A' (where K is the number of the separable operators) and hence 
is quite low. 

• The number of RPA states \v > is equal to the number of the relevant 1 ph configurations 
used in the calculations. In heavy nuclei, this number ranges the interval 10 3 -10 4 . Every RPA 
state \u > is characterized by the particular set of the values q u sk and p u sk which, following 
EH, self-consistently regulate relative contributions of different time-even and time-odd 
operators of the residual interaction to this state. 

• Eqs. EH, EH, EH), EH relate the basic SRPA values with the starting functional and 
input operators Q s k and P s k by a simple and physically transparent way. After choosing the 
initial operators Qsk, all other SRPA values are straightforwardly determined following the 
steps 


Qsk (0|[<5 s fc, Jf]|o) 


Y a 


sky 


^lsk,s'k' 


Psk ( 0 | [P s ki J s ] | 0 ) = 4 * X s k, ^sk,s'k' ‘ ( 82 ) 
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As is discussed in Sec. eh the proper choice of Q sk is crucial to achieve good convergence 
of the separable expansion © at a minimal number of separable operators. 

• SRPA restores the conservation laws (e.g. translational invariance) violated by the static 
mean field. Indeed, let’s assume a symmetry mode with the generator P sym . Then, to keep 
the conservation law [H, P sym ] = 0, we simply have to include P sym into the set of the input 
generators P sk together with its complement Q sym = i[H, P sym ]• 

• The basic SRPA operators can be expressed via the separable residual interaction (J2J): 

Xsk = -ifc e s p , P sk U tk = -i[VZ P , QskU (83) 


where the index ph means the 1 ph part of the commutator. It is seen that the time-odd 
operator P sk retains the time-even part of to build X sk . Vice versa, the commutator 
with the time-even operator Q sk keeps the time-odd part of V)|® p to build Y sk . Equations dHTH) 
hints also the relation between the SRPA operators and the true (not separable) residual 
two-body interaction. 

• Some of the SRPA values read as averaged commutators between time-odd and time-even 
operators. This allows to establish useful relations with other models. For example, m , 

(J22I) and m give 




-(o|[P,. t .,[c/,p st ]]|o), 

-(o|[4'*',Ky,e, t ]]|o). 


(84) 

(85) 


Similar double commutators (but with the full Hamiltonian instead of \Pg P ) correspond to 
ms and mi sum rules, respectively, and so represent the spring and inertia parameters Q 
in the basis of collective generators CV , and Psk- This demonstrates the connection of the 
SRPA with the sum rule approach 3S| and local RPA B- 

• The commutator form of the SRPA values allows to represent them through the matrix 
elements from the operators entering the commutators. Namely, the strength constants and 
responses gain the form (see Appendix IB 31) 


Ki,Kj> 0 

i < 0|[Q sfc ,R s ]|0 > = -4 y <ij\Q s k\0> S{< ij\B s \0 >} , (86) 

ij 

K u Kj> 0 

i < 0|[P S fc, -Rs]|0 > = -4 i ^ < ij\p ak \0 > 5P{< ij\B s \0 >} (87) 

ij 
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where B equals to X for the strength constants and to J (the density operator) for the 
responses. Since the involved matrix elements are already in our disposal, these forms 
considerably simplify the calculations. Besides, these forms are convenient for analysis. For 
example, they allow to determine the conditions of vanishing some response components, 
see Appendix IB 31 for more details. 

• In fact, SRPA is the first TDHF iteration with the initial wave function (HU). A single 
iteration is generally not enough to get the complete convergence of TDHF results. However, 
SRPA calculations demonstrate that high accuracy can be achieved even in this case if to 
ensure the optimal choice of the input operators Q s k and P s k and keep sufficient amount of 
the separable terms (see discussion in Sec. El In this case, the first iteration already gives 
quite accurate results. 

• SRPA equations are very general and, after simple modifications, can be applied to diverse 


systems (atomic nuclei, atomic clusters, etc.) described by densi 


functionals, see for the case of atomic clusters Refs. 


23J, 


y and current-dependent 


m 
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2.6]. Even Bose 


systems can be covered if one redefines the many-body wave function flU exhibiting the 
perturbation through the elementary excitations. In this case, the Slater determinant for lph 
excitations should be replaced by a perturbed many-body function in terms of elementary 
bosonic excitations. 

III. SPECIFICATION FOR SKYRME FUNCTIONAL 


A. Skyrme functional and its mean field Hamiltonian 

We use the Skyrme functional juj in the particular form HQQ 

E j dr in T Ffsk(/ts, T~s j &sijsi Js) T T~tc (.Ppj'j ddj C mi 


where 


B - —T 

'^kin 0 

lm 

n c = X [ df *Pp(r) I ~ 1 -?| Pp{?) - 7 ^Q f Mp)] f , 


( 88 ) 


3 2,3, 


I r — F' 


«sk = "X 2 - ^ - |p(Ap) + | £ p,(Ap.) 


4 7 r 

h 
2 


(89) 

(90) 

(91) 
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3 ck+2 


K 




+bi(pr - j 2 ) - b[ J^(p s r s - f s ) 

S 

-b A (p(V3) + a ■ (V x J)) -b' 4 J2 (ps(V$s) + <? s • (V x y s )) 

s 

are kinetic, Coulomb and Skyrme terms respectively. The densities and currents used in 
this functional are defined in the Appendix O Densities without the index s involve both 
neutrons and protons, e.g. p = p p + p n . Parameters b and a are fitted to describe ground 
state properties of atomic nuclei. The value E cm is the center mass correction. 

First functional derivatives read 


UJr) = 


SE 

bps(r) 


= bop{r) - b' oPs {r) 


(92) 


+ h 


a 


~P a+1 (r) - ^{ap a l {f)^p 2 Ar) + 2p a (r)p s (f)} 


3 ' w 3 

S' 

+ bir(r) - b[r s (r) - b 2 A? p(r) + b' 2 A? p s (r) 
— b A Vr- A(U) — b A Vf 5s(r) 

+ Kv f - (|) I/3 IP P (fO] 1/3 >, 


n 


7T 


BJr) = 


5E 


2m 


+ bip{f) -b\p 8 (f), 


5r s (r) 

SE 

W s (r) = —- = 6 4 Vr-p(r) + ^V^-p^r). 

<53 s (r) 

for time-even densities and currents and 
5E 


(93) 

(94) 


A,(r) = 


S.if) = 


bjsif) 

5E 

5a s (r) 


= —2b\j{f) + 2 b[j s (r) - b 4 (W P x a(r)) - 64 (Vf x a s (r)) , (95) 


= - 64 (V^ x J(r)) - 64 (Vf x j s (r)) 


(96) 


for time-odd ones. The last line in m is the Coulomb contribution. It includes only proton 
density (s — p). Since static ground-state time-odd densities are zero, the values a(r) and 
j[r) in (1951) and (1961) are reduced to the time-dependent density variations ©• We assume 
that the functional derivatives and the involved densities are have time dependent 

but, for the sake of simplicity, do not depict this. 

Using we get the time-even part of the mean held Hamiltonian (j8l) 

h s + (r) = U s {p , z) - VH s (p, z) V - iW s (p, z) • V x S (97) 

= U s (p, z ) + i VB s (p, z)V - iW s (p, z) -V x a. 
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If to implement the static densities and currents, one gets the ground-state Hamiltonian, i.e. 
h s + —> Iiq. The ground-state mean hied of even-even axial nuclei has no any contributions 
from the time-odd variables. 

The part of the mean held Hamiltonian (JHJ) from the time-odd densities reads 

ht[r) = - % -{A s {p,z), V} + )^{S s (p,z),a} (98) 

= • A s (p, z) - A s (p , z) ■ V) + S s (p, z ) ■ a. 

It is used only for derivation of the response Hamiltonian. Hence the involved current and 
spin densities are represented merely by their time-dependent variations. 


B. Time-even response 


For time-even densities, the second functional derivatives read 


8 2 E 

bU s (f) 

bp si (n)bp s {r) 

bps An) 


b 0 6„ 


(a(a- 
31 3 


a a V 

°S,S1°S,P 0 V 

O 7 

^S.Si^S.p I 

| r ■ 

5 2 E 

bU s {r) 

Sr 8l (fi)Sp s (r) 

br Sl {fi) 

8 2 E 

bU s {f) 

b% si (ri)bp s (r) 

^si(Fl) 

8 2 E 

bW s {r) 

5psi(ri)6%(r) 

bp si {fi) 


(99) 


s 2 


- 6 , 1 5 ss /]5(r - ri) , (100) 

i + b' 4 8 SSl ]X^b(r - rl) , (101) 

b 4 + b 4 S SSl ]W^(r - fi) . (102) 

The last two terms in m represent the exchange and direct Coulomb contributions. The 
pairing second functional derivative is not presented here. The pairing contribution to the 
response is considered in Sec. Ill Gl 

Following m, the operator Xf k (r) reads 

f X2 c 

= C; tl (r)p,(f) + B‘ tltl (r)f,(r) + W7 lh (f)S,(f) + S w x£‘ ,r) (103) 


x likSf) = 1 
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and its matrix element is 


< ij\X‘ lkl \6 > = J df{[/,V I (r>* i (r) + BJ lh (rX(f) + ^(r^r)} (104) 

+ i„, < > 

where p£-(r), rfj(r), and are transition densities (ID 3D - (ID 51) . Expressions for opera¬ 

tor A^ air) and matrix element < ij|X^ a "’|0 > are done in (linl and the Appendix 1(1 21 
respectively. The functions k (r), B* k (r), and k (r) read 

5 2 E 5 2 E 


= / dfi[ 


+ 


5p si (ri)5p s (r) 

5 2 E 


PX,siki(^l) + 


^x.sifei (d)] 


5r Sl (fi)<5p s (r) 


T x,siki (d) 


(105) 


53f si (fi)5p s (r) 

= U s Slkl (p,z) cos pd, 
r 5 2 e 

= / rfri = B s siki (p,z) cosp9, 


5p si (ri)5r s (r) 

-> f 5 2 E 

W° lkl {f) = / dr=i [—-+ 


<5 2 £ 


dp Sl {fi)5^ s (r) 

5 2 E 


ST Sl (r l )S^s s (r 


-'Tx,s 1 k 1 {f l) 


(106) 

(107) 


■Sxsifci (ni)] 


5S sl (fi)5Qy(r) 

= e p W^ lkl . p (p, z) cos pd + e z Wg lkl . z (p, z) cos p6 + e e W s s lkl . d (p, z) sin p0 

where the values U s sikl {p,z ), B s siki (p,z), W s sikl . p {p, z), Wf lkl . g (p, z), and W s s lkl . 0 (p, z) are 
given in the Appendix 1(121 

Contribution of the Coulomb integral to the matrix element m reads 


.,A, , 2 /' _ f, _ PX kl) p(Pl, Z l )COSp0 1 

<A X s 1 k 1 \‘]>Coui = d SUP Ss, P e j dr j dr i- ^ ^ , - Pijif) 


r — rp 


(108) 


= I dr Vh^iPi*) cos 


where U k ^ ul (p,z) is determined in (IF9D . It should be included to U* k (p, z), see 

In our case, the singularity problem for the Coulomb interaction cannot be treated by the 
familiar methods. These methods mainly deal with spherical systems Jj 34 | or axial systems 
in the static ground state 0 and, in any case, do not master the specific structure of the 
dynamical response in axial systems, given in (11081) . Hence we propose in Appendix IF1 the 
new prescription. 
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C. Time-odd response 


For time-odd densities, the second functional derivatives read 

= 2 [-&! + &i<S a , ai ]5(fi - r) , 

= ~[h + 64 ^, 81 ]x <5(ri - r) , 

,, ^ , = -[ & 4 + x «5(fi - r) . 

0JsA r i)0(Ts{r) 6j Sl [r 1 ) 

Following (jTTflt . the operator (r) has the form 

5 2 £ 


§ 2 E 

5A s (f) 

Sj si (ri)Sj s (r) 

fijs 1 (fl) 

5 2 E 

SA s (r) 

5a Sl (fi)Sj s (r) 

Sd si (f 1 ) 

5 2 E 

<s£(f) 


(109) 

( 110 ) 
(111) 


Y s \kA?) = * 


i f dr 1 

aa-L 


1[ 6J«AA)6J?(f) 

= A s sikl (f)j' s (r) + S* aikl (r)a a (r) . 
and its matrix element is 


j,J<0|[Q, lfcl ,J«T]|0>(f 1 )J“(r|] 


( 112 ) 


< 'ilhUJO >= / +s; itli ,(f)sy i# (f)] (113) 


where 




<5 2 £ 


'5j S i(ri)(5i s (r) 


jV,Slfcl + 


<5 2 £ 


5<? si (ri)hj s (r) 


■Sy,siA:i] 


2[ b\ “b &i^s,si]jV,sifci [^4 4~ ^4 ( 5 s,si](V X Sy^j^) 

= epA a Slkl . p (p, z) cospd + e z A s Siki . z {p, 2 ) cos pQ + e e A s Siki . e (p, z) sin p6 , (114) 

6 2 E 


S s slkl (rA = dfi 


fijsi (ri)5a s (r) 


jY,siki 


[&4 d - ^4^s,si](^ X jy.sifei) 

= e p 5| lfci;p (p, 2 ) sin p9 + e z S s slkl . z (p , 2 ) sin/r0 + e^ lfel . e (p, 2 ) cos^0 , (115) 


and 


rP 


(116) 


(V X Sy jS fe) = e p [-sy iSfe ; Z (p, 2 ) - d z S Y ,skAPi z)] cos p6 

4“ C z [— dp{^pSY t sk;o') ~'5y,sfc;p] COS pd -(- eg[^Sy" s fejp <9p<Sy,sfc;z] sill [id , 

(V Xjy jSfc ) = e p [--j Y>sk . z (p, z) -a 2 jy )Sfe;0 (p,2)]sin^6> (117) 

+ e z \^d p (pj Y ,sk-,e) 4 - ^jy, s fc;p] sin/r# + e e [d z j YjSk . p ~ d p j Y>sk . z \ cos pd . 
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Expressions for A^ sk . p (p, z), A^ sk . z (p,z), A^ ak . 0 {p,z), S^ sk . p (p,z), S^ sk . z (p,z), and 
Sy s k-e{Pi z ) are given in the Appendix 1(1 31 


IV. CHOICE OF INITIAL OPERATORS 


As was discussed in Sec. ITITI all SRPA operations start from initial (generating) operators 
Q s h, see the sequence of the model steps in ( 1821 ) . The SRPA formalism itself does not provide 
these operators. At the same time, their proper choice is crucial to get good convergence 
of the separable expansion 0 with a minimal number of separable operators. The choice 
should be simple and universal in the sense that it can be applied equally well to all modes 
and excitation channels. 

We propose the choice inspired by physical arguments. The main idea is that the gener¬ 
ating operators should explore different spatial regions of the nucleus, the surface as well as 
the interior. The leading scaling generator should have the form of the applied external field 
in the long-wave approximation, which is most sensitive to the surface of the system. Since 
nuclear collective motion dominates in the surface region, already this generator should pro¬ 
vide a good description. Next generators should be localized more in the interior to describe 
an interplay of surface and volume vibrations. For AA-excitations in spherical nuclei, the 
best set of the generators was found to be [nj] 

Qk(f) = Rk(r)(Y Xp (n) + Y* p (n)) , P k = i[H,Q k ] (118) 


with 

r A , k — 1 

j\(q k \r), k = 2,3,4 

q\ = a k ttK a 2 = 0.6 , a 3 = 0.9 , a 4 = 1.2 

Adiff 

where Rdifr is the diffraction radius of the actual nucleus and za is the first root in j\{z\) = 0. 
The separable operators X k and X k with k — 1 are mainly localized at the nuclear surface 
while the operators with k > 1 are localized more and more in the interior. This simple set 
seems to be a best compromise for the description of nuclear giant resonances in light and 
heavy nuclei. 

In deformed nuclei, we exploit even simpler set with 

R k (r ) = r A+2(fc - 1} . (120) 
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Similar to the previous set, the separable operators X k and X k with k — 1 are mainly 
localized at the nuclear surface while the next ones are localized more and more in the 
interior. We expect that already the first generators with k — 1 and 2 are quite enough for 
description of giant resonances. 

In the case of magnetic modes, the initial generators should be the time-odd. Following 
the same logic, the leading (k=l) operator should coincide with the transition operator for 
the external magnetic field in the long-wave approximation. The next generators can be 
produced following the prescription (11201) for the radial parts. The time-even conjugates of 
the time-odd generators can be obtained as Q k = i[H, Pk] p h- 

V. CONCLUSIONS 

A general procedure for self-consistent factorization of the residual nuclear interaction 
is proposed for arbitrary density- and current-dependent functionals. The separable RPA 
(SRPA) constructed in the framework of this approach can dramatically simplify the cal¬ 
culations while keeping high accuracy of numerical results. The economical effect of SRPA 
is especially actual for deformed nuclei. In the present contribution, SRPA is specified for 
description of axial nuclei with Skyrme forces. 

SRPA can be used for description of EX (and MX) response in both spherical and de¬ 
formed nuclei. The approach can also serve for getting the basis of one-phonon RPA states 
for further description of anharmonic corrections, vibrational admixtures in the low-energy 
states in odd and odd-odd nuclei, etc. One of the most promising lines of future studies is 
dynamics of exotic nuclei obtained in radioactive beams. 
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25 


APPENDIX A: CYLINDRICAL COORDINATES 

1. Expressions for cylindrical coordinates 

Cylindrical coordinates are defined as 

x = pcos'd e x = e p cosd — e^sim) e p = e x cosd + e y sin%) 

y = psind e y = e p siri'd + e$cosd e$ = —e x sim) + e y cosd 
z = z e z = e z e z = e z 

The gradient operator reads 

V = e p V p + e^V^ + e 2 V z 


where 


1 


iL z 


V p = d p , Wo = -do = —, V 2 = d ; 


P P 

and L z = — ido is the third component of the orbital momentum operator. 
The Laplacian, divergence, and curl read 


A — V • V — d p H— d p H —-dfl + d^, 

P P 

divA = V • A = —d p (ypAp) H —doAo + d z A z , 
rot A = \^d#A z - d z Ao\e p + [d z A p - d p A z \eo 
J r[—dp{pAo) — —doA p ]e z . 


The vector of Pauli matrices is 


(j c p a p T eo (Jo T e 2 c z , 


where 
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2. Single-particle wave function in cylindrical coordinates 


The single-particle particle wave function and its time reversal are expressed in the cylin¬ 
drical coordinates as spinors 


^i(r) 


R^\p, 


(Al) 


%(r) = f &i(r) 


r\ + \p, z)e~ im i +) ^ ) 


(A2) 


where K t is the projection of the complete single-particle moment onto symmetry z-axis of 
the axial nucleus. 

In short notations covering both the state i and time-inverse i states, the spinors read 


a) (p,z)e 


(< t ) q 

2771 - '77 




R\ (p,z)e zmi 
-aR { ~ a \p , z)e- imi ~ a) ^ 


for ordinary state i 
for time reversal state 


where 


(<t) 7V- ^ 

m) — K -cr, 

1 2 


(—cr) (cr) . 

m- = m\ + u, 


(cr) (— or) 

rn- = —m) . 

I r, 


(A3) 


(A4) 


APPENDIX B: USEFUL RELATIONS 

1. Hermitian and time-conjugate properties 

All the operators used in the model are hermitian (A = A^) and have the definite time 
density 

A = TAT- 1 = 7^A, 7^ = ±1 

where T is the time-inversion operator: T\i >= |i >, T\i >= — \i >. 

There are useful transmutation relation for the hermitian operators 

< j\A*\i >=< j\A\i >=< i\A\j >* (Bl) 

and relations for time-inverse states and operators 

< j\A\i > = <j\T- 1 TAT- 1 T\i>=<Tj\TAT- 1 T\i>’ t =<j\A\i>’*, 

< j\A\i > = - < j\A\i >* . (B2) 
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In axial even-even nuclei, the densities are built from the pairs 


< i\A\i > + < i\A\i >=< i\A\i > +7^ < i\A\i >*=< i\A\i > (1 + 7^) (B 3 ) 

and hence vanish for time-odd operators. 

2. Connections for operators with definite time-parity 

Operators with the definite time-parity have the useful property 

< 0|[i,R]|0 >=< 0|AB|0 > (1 -7t7t) ( B4 ) 

i.e. the average value of the commutator survives only if operators A and B are of the 
opposite time-parity. Indeed, 

< 0|[i,R]|0 >=< 0|AB|0 > - < 0|M|0 > 

= < 0\AB\0 > - < 0 \T- 1 TBT~ 1 tAt- 1 T \0 >=< 0\AB\0 > -7 £7$ < Q\T- x bAt\Q > 

= < 0|AB|0 > -7^71 < T0\bA\T0 >*=< 0|AB|0 > ~^t!t < 0|si|0 >* 

= < 0 |iS |0 > -7^7r < 0 |( 5 i ) + |0 >=< 0 \AB \0 > (1 - 7^77) . (B 5 ) 

3. Representation through matrix elements 

Responses m-m and inverse strength matrices m-m read as the averaged com¬ 
mutators 

( 0 |[R,R]| 0 ) with 7 r = — 7 t- (B6) 

Calculation of these values can be greatly simplified if to express them through the matrix 
elements of the operators A and R, because in practice these matrix elements are already 
in our disposal. In this case, the commutator reads 

< 0|[i4 a ,5 s ]|0 > = y^{< 0\A s \ij >< ij\B s \0 > - < 0\B s \ij >< ij\A s \0 >} (B 7 ) 

ij 

Ki,Kj> 0 

= 2 ^ {< ij\A s \0 >*< ij\B s \0 > - < ij\A s \0 >< ij\B s \Q >*}. 

ij 

The sum runs all the states, both ordinary and time reversed. Using the properties m- 
(IB2I) . the matrix elements < ij\A s \{) > and < ij \ A s \0 > are reduced to < ij \ A s \0 > and 
< ij\A s \Q >, respectively. Hence the coefficient 2. 


In the SRPA, the operator A is associated with operators Q or P, which have real or 
imagine matrix elements, respectively. Then 

Ki,Kj> o 

< 0 \[Q sk ,B s \\0 > = 2 22 < ij\Q s k\0 > (< ij\B s \0 > - < ij\B s \0 >*) 

ij 

Ki,Kj> 0 

= 4 i 22 <ij\Qsk\0> 3{< ij\B a \0 >} , (B8) 

ij 

Ki,Kj> 0 

< 0\[P sk ,B s ]\0 > = -2 22 < ij\P s k\0 > (< ij\B s \0 > + < ij\B s \0 >*) 

ij 

K i} Kj> 0 

= -4 y t <ij\P sk \0> 3?{< ij\B s \0 >} (B9) 

ij 

where A{...} and mean the imagine and real parts of the values in the parenthesis. 

Then the strength matrices read 


^s'k' ,sk 


^Is'k' ,sk 


Ki,Kj> 0 

-i<0\lP s , k ',X s s ' k \\0>=4i <ij\p s , k ,\0> $t{< ij\X s s ' k ]\0 >} 

ij 

Ki,Kj> 0 

-4 22 < ij\P s 'k'\0 > < ij |X|']|0 > , (BIO) 

ij 

Ki,Kj> 0 

-i(o|[g fc ;,n]|o>=4 22 < i J\QMo> 3f{<u|y/ fc ']|o>} 

ij 

Ki,Kj> 0 

4 < i j\Qs'k'\0> <ij\Y s H\0> (Bll) 

ij 


where the overline matrix elements 


<ij\P s 'k'\0> = -i <ij\P s 'k'\0>, < ij\X°' k ]\0 > =< ij\X? k ]\0 >, (B12) 

<ij\Q s 'k'\0> = <ij\Q s 'k'\0>, < ij\Y s s p \0 > = -i < ij\Y ? k ]\0 > (B13) 

are real. It is seen that both strength matrices (IB 101) and (IB 111) are real as well. 

The case of responses is more involved in the sense that matrix elements of the second 
operator in the commutator are transition densities (ID II) - (ID 51) which are generally complex. 
The responses entering X and Y operators read 

Ki,Kj> o 

n\i sk = i < o|[Afc,i“]|o >= ~4i 22 < y|Afc|o > K{< ij\j?\o >} 

ij 
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Ki ,Kj> 0_ 

= 4 < ij\P s k\0 > 3?{< ij\Jg\0 >}, (B14) 

ij 

Ki,Kj> 0 

= i<0|[Q sfe ,J s a ]|0>=-4 ^ <ij\Qsk\0> %{< ij\J a \0 >} 

ij 

Ki,Kj >0_ 

= -4 < ij\Qsk\0 > %{< ij\J?\0 >} (B15) 

ij 

where < ij |,/“ 10 > are transition densities. It is seen that all the responses are real. 

Following m-m, the transition densities (or their components) for p = K, ± Kj are 
proportional to e~ = cos pi) — i sin pd or = i cos pd + sin pd. Hence the response 

obviously vanishes at p = 0, if its £y{...} or part delivers sin /id. 


APPENDIX C: DENSITIES AND CURRENTS FOR SKYRME FUNCTIONAL 


In Skyrme forces, the complete set of the densities involves the ordinary density, kinetic- 
energy density, spin-orbital density, current density and spin density: 


Ps(r,t) = ^2(p* h {r,t)(p h {r,t), 
hes 
occ 

Ts(r,t) = ^2Vip* h (r,t)-V(p h (r,t), 

hes 

occ 

(r,t) = ~i^2<p* h (r,t)(V xa)(p h (r,t), 




hes 


3s(r,t) = ^J2\Ph(r,t)Vp h (r,t)-Vv* h (r,t)ip h (r,t) 


hes 


&s (r) = 


^2<Ph(r,t)#<p h (r,t), 


TpT- 1 = p 


Trf- 1 = t 
TSf- 1 = 3 

TjT - 1 = -j 
Taf ' 1 = -a 


hes 


where the sum runs over the occupied (hole) single-particle states h. The associated hermi- 
tian operators are 

N a 


Ps(r) = J^Sifi-r), 
i =1 
Ns 

^V<5(r)-r)V, 


t, r = 


i= 1 
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3 s(r) = -r)Vxa, 

i =1 

3 sir) = 

i— 1 
N s 

& s {r) = ^<5(r; - r)a . 

Z— 1 

where a is the Pauli matrix, N s is number of protons or neutrons in the nucleus. 


APPENDIX D: TRANSITION DENSITIES 


1. Explicit expressions 


For p = |Ki — Kj | with K t , Kj > 0, the transition densities read 


< 

Ws 

|0 

> 

(7) = 


■ P iPM 

) + 

• Via 

zip, 

Z ) + e# 

■fij-fp, z)\e 

-i(I<i-KM 

, 

(Dl) 

< 

ij\s s 

|0 

> 

(f) = 

= fp'fr, 

P ip,z) 

+ e z • 

4; * 

(P,z 

) + e# ■ 

ifrfP, 

))* 

-i(Ki-KM 

, 

(D2) 

< 

ij\Ps 

|0 

> 

(f) = 

= Pij(p,z) 

e~« Ki 

-KM 







(D3) 

< 

ij\r s 

|0 

> 

(f) = 

TT 

r* 

e-idu- 

-KM 







(D4) 

< 

ij\% s 

|0 

> 

(r) = 

= [e P -% 

■pip,Z 

) + 


ziP, 

z) + e,? 

■ fp, 

z)] 

e -i(Ki-I<M 

(D5) 


where the components with low indices ij are real. The case of the time-inverse state j is 
straightforwardly obtained by Kj —> —Kj. As is shown below, the components accompanied 
by the imagine unit (jfj. p ,jfj. z , sf r g , and 3f|-. e ) demonstrate the specific features. 

Using wave functions from the Appendix IA 21 and density operators from the Appendix 
□ one gets expressions for the real components of the transition densities m-m and 
for their transmutation properties: 


itj- P ip, z ) 

fir, Z {P, *) 

fij-fp, z) 


\K ] e 

i 

^ s 

1 

^ s 

_1 

cr=+,— 


\KY. 



(T= + ,— 




r?U) SU) 

n i n j / — (cr) . ~{aY 

-— (m) + m) ' 

p 3 


(D6) 

(D7) 

(D8) 


fji:n(P , Z ) = fj-piPiZ), jji-fp, Z) — z (p, Z ) 


3ji;*(P’ Z ) = -firfP’ Z )■ 
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sij-p(p,z) : 

Sij;z(p, Z) 
S ij; i?(P> Z ) 


E [tf 


«<"> V I ¥&- < ’ ) R { ° ) 

L J 


4' 1 E 


E 


cr 


a 


p(°9 p( CT ) 

ji, ttj 


h(~ a ) p( CT ) 

K i K j 


(D9) 

(DIO) 

(Dll) 


s ji-,p(p, z ) = ~ N *j: M z ), S jr,z(P, z ) = ~S S ij- z {pi Z), S s ji; #(p, Z )=S s ij .#(p i Z). 


p s ij(p , z ) = u { 


- 7/(+) 


*7 z_^ 




(D12) 


pUp, z ) = pUp, z )• 


r?.(p,«) =z4 h) 5] (d p R^)(d p R^) + (d z R^)(d z R^) + 


rDE m i C>D) dD) 


P^ 


TD\ 

- Jh^ rij 


(D13) 


%- tP (p, z ) = 


z ) = ^(p, -)• 

V E 


a 


E 


R^\d z R ( r a) ) + Rf\d z R ( - a) ) 

p(°0 p( ff ) 

, / ~M . ~ (ct)x -“■i n j 

+ (m- +m) ; )- 

P 

a(R { r c \d p Rf ) ) + R { r a \d p R ( ^ 


p( CT ) p( —<T ) 
/ ~ (—cr) ~ (<r)\-*h 

-(m\ ' + m) ) --— 

P 


(D14) 


(D15) 


% ; *(p,*) = 7^ ^ [a(^ ) (^ ) ) + J Rj.-- ) (^- CT) )) 




(D16) 


%% 0 (p,z) = 


ji; 


— 

*j; p 


(P, *), 


cv 


(P, 


%;*(P. *)» Z) 


~%-;AP, Z )' 


2. Useful features 

Some components of the vector responses become zero at p — 0. These cases can be 
revealed by inspecting the specific combinations of the transition densities involved into 
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responses. Following (IB 141) - (IB 1 5(1 . these combinations are 


< ij\Jg\0 > + < ij\Jf\Q >*= 2 K{< ij IJ“|0 >} 
for time-even densities 

< ij |J“|0 > - < >*= 2i 3{< ij\ J S Q |0 >} 

and for time-odd ones. 

Specifically, these combinations read 

< ij I/JO > - < ij\j s \0 >* 

= 2 *[(e P • jtj. P {p, z) + e z - jf j; z (p, z )) cos/id - e# ■ j? j; #(p, z) sin/id], 

< r/|s s |0 > - < ij\s s \0 >* 

= 2 i[-{e p ■ sfj. p {p, z) + e z - sb z (p, z)) sin /id + e* • sb #{p, z) cos//??]. 

< ij\ps\Q > + < ij\ps\Q >*= 2 p s ij{p,z) cos pi), 

< ij\f s \0 > + < ij |f.s|0 >*= 2r-j(p,z) cos/I'd, 

< ij|3 s |0 > + < ij |3 S |0 >* 

= 2 [(e P • 3b p O, -) + e 2 • 3b *(p, z)) cos /id - • 3b „(p, 2 ) sin /id] . 

For /i = 0, the combinations with sin/id vanish 

(< r/|3 s |0 > + < ij|3 s |0 >*) e = 0, (< ij\j s \0 > - < ij\j s \0 >*)e = 0, (D17) 
(< ij\s a \0 > - < ij\s s \0 >*) p = 0, (< ij\s s \0 > - < ij\s s \0 >*) z = 0 (D18) 

and hence the corresponding response components 


d — 0 ~^ JY.sk — S Y,sk — S Y,sk ~ ^ Y,sk ~ 0- (^19) 

The curl in cylindrical coordinates is 

rotA = VxA = \^d$A z - d z A# ]e p + [d z A p - d p A z \e# + \^-d p {pA^) - -d#A p ]e z . (D20) 

Then, taking into account (inuD-dnuD and zero value of the derivatives dg from all the 
transition densities at /i = 0, one gets 


33 










These relations are used in derivation of the first and second functional derivatives 
of the Skyrrne functional (terms with 64 and b' 4 ). 


APPENDIX E: RESPONSES 

The responses have the general form (piiD - dHThl) . The explicit expressions read 

Jy,sk{r) = e p j^ sk (r) + e z j^ ak {r) + e e j^ sk (r) 

= e pjy )Sk (p, z) cos pd + e z jy sk (p, z) cos p9 + e g j^ sk (p, z) sin p9 , (El) 
s Y ,sk(r) = ^P s Y, s k(^ + ^Sy, sk (r) + e g s e YtSk (r) 

= Cp s Y,sk(P, z ) s[n P 9 + ZzSy tSk (p, z) sin p9 + egSy sk (p , z) cos p6 , 

Px,sk{r) = Px,sk(p, z) cos pd 
Tx,sk{r) = Tx,skip, z) cos pd 

3 A >fc(r) = e p 3 p x>fc (f) + e z ^ z XiSk [f) + e e ^ e x , sk {f) 

= ep^x,sk(Pi z ) cos /^ + e*Sx )Sfc (p, z) cosp0 + e e % d X sk (p, z ) sinp0 , (E2) 

The components of time-odd responses have the form 

Ki,Kj> 0 

j Y ,skif) = 1 < tt\[Q S k,js]\o >= -4 ^ < *j|<? S fc|6 > 3{< f.7(if]|o >} 

ij 

Ki ,Kj >0_ 

= -4 J < U'lQafclO > jtj. P {p, z) cos pd = jy jSk {p, z) cos p9, (E3) 

ij 

Ki,Kj> 0 

= * < Q\[Qsk,Js]\Q >= -4 ^ <ij\Qsk\0> 3f{<y|5s]|0>} 

ij 


34 






Ki,Kj> 0 

= -4 Y < ij\Qak\0 > jij. z iP, Z ) COS P° = JYskiP, Z ) COS pO, (E4) 
ij 

Ki,Kj> 0 

3Y,sh{r) = i < >= “ 4 Y <ij\Qsk\0> %{<ij\j 9 s ]\0>} 

ij 

Ki,Kj>0 _ 

= 4 Y <ij\Qsk\0> jij., eip, z ) sin pO = jf >fc (p, z) sin pO r (E5) 

ij 

Ki,Kj> 0 

s Y,sk(r) = 1 < 0|[Q»fc»s£]|6 >= -4 Y < ij\Qsk\0 > ^{< ij\s p s ]\0 >} 

ij 

Ki,Kj> 0 

= 4 Y < ul&fc|6 > 4; p(p> *) sin / i0 = s y, s fc(P z ) sin pO, (E6) 

ij 

Ki,Kj> 0 

4, s fe(^) = * < 0|[Q»fc,s*]|0 >= -4 ^ < ij\Qak\0 > ${< u|s*]|0 >} 

ij 

Ki,Kj> 0 _ 

= 4 < Vl4 fc |0 > 4; *(P> *) sin /^ = s Y,skiP, z ) sin pO, (E7) 

ij 

Ki,I\j>0 

SY, ak (r) = i < 0|[Q*fc,sf]|0 >= 2 i Y < ij\Qsk\0 > 2i ${< u|s®]|0 >} 

ij 

Ki,Kj> 0 

= -4 Y <i j\Qsk\0> sfj. g (p,z) cos p6 = s e Ysk (p, z) cos pO . (E8) 

ij 

The time-even responses and their components are 

Ki,Kj> 0 

Px,skir) = i < 6|[P S fc,Pa]|0 >= ~4i Y < y|[Afc|0 > 3?{< ij|p s ]|0 >} 

ij 

Ki,Kj> 0 _ 

= 4 < ij|P s fc|0 > p^(p, z) cos p6 = p x ,skip, z ) cos pd, (E9) 

ij 

Ki,Kj> 0 

Tx, s jt(0 = * < 0|[Afe,T s ]|0 >= -4 ^ < ij\[P sk \0 > K{< ij\T s ]\0 >} 

ij 

Ki,Kj> 0 _ 

= 4 Y < ul Afc|0 > r*(p, z) cospd = Tx,skip, z) cos pd, (E10) 

ij 

Ki,Kj> 0 

3y iSfe (r) = i < 0|[P sfe ,^]|0 >= —4i < y|[A fc |0 > K{< y|3£]|0 >} 

ij 

Ki,Kj> 0 _ 

= 4 < ij I Afc|0 > P (p, z) cospd = c y x ,skip, z ) cos pd, (Ell) 

ij 
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Ki,Kj> 0 


3 Z 


W*) = ' < 0|[-P.«=, 3J]|0 >= -4* £ < 0I[A*|8 > *{< ij|a;]|0 >} 




Ki,Kj> 0 


4 < u|Afe|0 > %. z (p, z) cos p6 = $Sx >ak (p, Z ) cospO, (E12) 


*7 


Cxd 

x) 


e y Jk (f) = i < 0|[p, t ,3f]|0 >= -4i V < ij|[P.»|0 > S{< V|9j]|0 >} 


Ki,Kj> 0 


-4 5^ < ^1^1° > %•; fl(P» *) sin l J ° = %x,sk(p, z ) sin p6. (E13) 

u 


The overlinc matrix elements are defined in (lFT2l) - (lBm The pairing factors m are 
included into the matrix elements and transition densities. The values are real. 

Time-even responses have diagonal contributions i — j while time-odd ones not. 

For p = 0, the responses fulfill (ID 191) and 



It is easy to see that the properties (nnsi) and (IE 141) of the responses fully repeat 
the properties (lDT7l) - (lDT8l) and (ID21I) of the combinations of the transition densities, 
entering the responses. 


APPENDIX F: COULOMB CONTRIBUTION 

The contribution of the Coulomb integral to the SRPA response is 

U£ oul (r) = e 2 [ dr 


rCoul / —i 2 / j-./'Vj.i 1 '/''-:')««/<«' 




(FI) 


The typical trouble is connected with the singularity at the point f = r'. Because of the 
angular dependence cos pd' in the nominator of (IF 1 1) . this trouble cannot be circumvented by 

n nn 

common methods. For example, the FTT solver method |24| and procedures (23,do not 
assume any angular dependence of this kind. So, we should develop our own prescription. 
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First, we bypass the logarithmic singularity by using the Vauterin identity [35l i36] 

A I- Wl 2 

A jf/\r — r = ———r. 


r — r 


(F2) 


Then, after integration by parts, the integral CD is reduced to 


Uf oul (r) = 


— J dr \r-r'\A^ p x , P k(p',z') cos pd' 
j J p'dp'dz'[d 2 p , + jd r/ - -^r 2 + d 2 z ,]p x , P k{p ', z') 

/* 2tt 

/ dd'\r — f '\cos pd'. 

Jo 


(F3) 


For p = 0, the angular integral in (El is recast to 

f dd'\r — r\cos pd'= A^/d(j)^z)E(—^-—) 


(F4) 


'd(p,z) 

where d(p, z) = (p + p’) 2 + (z — z 1 ) 2 and E( t p p p ^ ) is the elliptic integral of the second order. 
This integral can be approximated by a standard polynomial formula, ffowever, in the 
general case p > 0 we cannot get this result. So, we should develop another procedure. 

First of all, we should take into account that all the terms in <m have the common 
angular dependence cos p6. This should be the case for the Coulomb term as well. To prove 
this, it is convenient to rewrite the angular integral from m in the form 


r*27r 


I(r;p',z')= / d6'\r — r '\ cospd 1 = Ii cos/id — J 2 sin p6 

Jo 


(F5) 


where 


/»27T /»27T 

Ji(p, z; p', z') = / dd'\r —r'\cosp(d'— d), I 2 (p, z; //, z') = / dd'\f — r \ sin p{d' — d). 

Jo Jo 

(F6) 

One may show that I 2 = 0 and I\ does not depend on the angle d. ffence ESD is reduced to 


I — h cos pd 

thus delivering the desirable angular dependence in IF ID . 

Altogether, the Coulomb contribution Q to the response is recast to the form 


(F7) 


= cos/,0 


Coul ( 


(F8) 
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with 


1 


Uk OUl (p, z) = e 2 J p'dp'dz' h(p, z ; p', z') ■ [d‘ 2 , + ^dp - - + d 2 z ,\ p x , P k(p', z') 

= e 2 J p'dp'dz' j ddi[(p + p') 2 + {z — z') 2 — App' cos 2 di} 1 ^ 2 cos2/idi 

■[dp + ^ d P' - ^2 + d-T P fc(p', *') • 

The matrix element m reads as 


(F9) 


r2ir 


< > Coui~ e 2 / pdpdz Ug oul (p,z)p s ij (p,z) • / dO cos p9e Hl ° . 


*r 


(F10) 


APPENDIX G: MATRIX ELEMENTS 

This Appendix represents the matrix elements of the operators Psk, X s s ' k and Y s s k . The 
pairing factors are supposed to be included into the matrix elements. 


1. Matrix elements of operator Psk 

The matrix element of operator P s k reads 

< ij\P 9 k\0 > = 2i £ij < ij\Qak\Q > ~ < ij\Yfk\0 > (Gl) 

= i{( 2 e,j< ij\Q s k\0 > — < k/’lkj^lO >) (S^Ki-Kj + d-^Ki-Kj) 

+ < u|^7fe 2) |0 > (S^Ki-Kj - S^Ki-Kj)} 

for p = | Ki — Kj | and 

< ij\Psk\0 > = Zidj < ij\Qsk\0 > ~ < ij\Yf k \0 > (G 2 ) 

= i{26ij< ij\Q s k\0 > - < ij\Ysk l) \0 > + < u|^f |0 >}S^ Ki+Kj 

for // = Ki + Kj. Here, Ki and Kj are projections of the complete single-particle moment 
onto symmetry z-axis of the axial nucleus. Expressions for < ij\ Y* k ^ 1 0 > and < ij\ Y* k 2> |0 > 
are given in the section ITTdl The combinations of Kronecker symbols follow from the angular 
integrals. Because of the specific combinations in m, all the terms of this matrix element 
have the same permutation properties for i j. 
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2. Matrix elements of operator X 


sk 


The matrix element of operator X*' k reads 


< Ul^Cfcl 0 > - < + d^Kj-Ki) 


+ < y|^ s S fc (2) |0 >(^,K i -K j - 5^ Kj -Ki), 


(G3) 


for n = \Ki — Kj | and 


< ij\Xs k \0 >= {< ij\X s s ' k {1) \0 


ij\X s J 2) \0 >}5^ Ki +Kj 


(G4) 


for /i = A'j + Aj. 
Here 


<yK |0> = Try pdpdz [U s sk (p,z)p s ij (p,z) + H s s fc (p, 0 )r^(p,^) 
+ A Wf k . g (p, Z)S<. Z (P, Z)] 

+ 6,s< tj\x:i {pair) \o > 


<ij\X% 2) \0> = 7r / pdpdz W&. 0 (p, z)Q( j;e (p, z) , 


<y|^- ) |0> = S^oSijGsXx, ks(uf - «*) ■ 


(G5) 


(G6) 

(G7) 


The values entering (1(471) are defined in Sec. mm 
Further 


Usk(px) = {bo~b' 0 S s 


+ [— b 2 + b' 2 5s,s'] [d 2 p + — ^7 + <%] + -7T---p“(Pj 

P p- 6 


- u 


, M a - !) «- 


p a 2 {px)J2 p2 s”(p’ z '> 


(G8) 


+ (p, z)\ps(p , z) + Ps'(p, -)] + -p“(p, *)<*-, *'] 

- S s >, p 8 sp ^-{—) 1/3 Pp 2/3 (p, z)} Px,,fc(p, z) + 8 S > p 5 sp ul Coul \p, z) 

3 7T 

+ [6i - &i<5 a , a /] Tx,sk(p,z) 

- [64 + 64^,-'] [^p(p^x, s fc(p> -)) + dz$fx, 8 k(p, z )) + ^x, 8 k(p > *)] > 


B s sk{pi z ) = [bi - b[S s ^}px,sk(p^ z ) ^ 


(G9) 
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(G10) 

(Gil) 

(G12) 


W? k ., p (p,z) = [64 + b' 4 8 s , S ']d p p x ,sk(P, z ) > 

W?*; *(/>, Z ) = [&4 + K S s,s']d z pX,sk(Pi Z ) » 

= -[&4 + &4^,s']^PA>fc(p,z) 

The overline matrix elements are real. The combinations of Kronecker symbols follow from 
the angular integrals. Expressions for f/^ Cou1 ^ (p ; z) in (IG8D is done in (IF9I) . The matrix 
elements (EH and m vanish at p = 0. Instead, the pairing matrix element O exists 
only at /i = 0. 

3. Matrix elements of operator 


The matrix element of operator Yf k reads 

< U’l^7fc|0 > = i{< ij\Y s s k {1) \0 >(6^-^ + S^Kj-Ki) 

~ < ij\Y s l {2) \U >(^,K i -K j - (G13) 

for /i = \K t — Kj | and 

< v\Yi\0 >= i(<WY¥> - < WSKi >K,k,+k, (G14) 

for p = Ki + Kj. The over line matrix elements are real. 

Further 

< yl>l (1, |6 > 


< »\y;Y\o > 


^k;p(P> *) 

2[ ^1 H - jY,sk;p{,Pi 

(GIT) 


[&4 “I - ^4^s,s'] dzSy^k-O^p^ %)\ 7 


A s s ' h ,z(P, Z ) 

2[ &1 + ^i^s,s']jX,sA;;2: 

(G18) 


= 7 t J pdpdz (G15) 

• Wsk-AP' Z )dij;p(P’ Z ) + z )jiy,z(P’ Z ) + Z )]> 

= 7 t J pdpdz (G16) 

• [^fe;o(P. z )ji-,e(Pi z ) + *)]• 
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(G19) 


[^4 “1“ ^4 ^s,s'\\^pi^P^Y,sk\d) ^Y,sk\p[ 

A s s k , e (p,z) = 2[-6i + V 1 5 a ,s'\jY,8k-fi 

[&4 + ^4^s,s'] [9z^Y,sk-,p dpSY,sk-,z\ i 
S S Lo(pi Z ) = h + b 4 $s,s'][-jY,sk;z{P, z) + d z j Y , a k-fi{P, z)\ , 

P 

sk;z(Pi z) [^4 T b^S s ,s'}[~9p{pjY,sk-,d) 4" ~jV,sfc;p] j 

sk;oiPi z ) [^4 “I" b ^s,s'] [^zjY,sk;p 9pjY,sk;z\ ■ 

For n — 0, one puts 

KM Z ) = S*;p(P> z) = S S sl,M z) = 0 
since these terms adjoin sin jid. The remaining terms are reduced to 

A S s kl p(P, Z) = 2[-&i + b'^s^j^sk-Ap, Z) - [64 + b' 4 5s,s'] d z S Y,skAP, z) , 
^~sk',z(.Pi z) 2[ bi T b^8s,s'\jY,sk-,z [^4 T ^4^s,s'] ~^pi,P^Y,sk',o') 1 

S s k;d(Pi z) = — [ b 4 T b 4 5s,s']\PzjY,sk;p ~ dpjY,sk;z] ■ 


(G20) 

(G21) 

(G22) 

(G23) 

(G24) 

(G25) 

(G26) 
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